function Jacobian(A,y)
len = length(A);
B = zeros(len,len+1);
for i = 1:len
    for j = 1:len
        if i == j
            B(i,j) = 0;
        else
            B(i,j) = -A(i,j)/A(i,i);
        end
    end
    B(i,len+1) = y(i)/A(i,i);
end
X = zeros(100,len);
for i = 2:100
    for j = 1:len
        X(i,j) = sum(X(i-1,:).*B(j,1:len)) + B(j,len+1);
    end
end
disp(X)
end